NAS A Technical Memorandum 106585 


/SS-3? 

■ygaay 


Comparison of Life Theories for 
Rolling-Element Bearings 


Erwin V. Zaretsky 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

Joseph V. Poplawski and Steven M. Peters 
J.V. Poplawski & Associates 
Bethlehem, Pennsylvania 


Prepared for the Annual Meeting 

sponsored by the Society of Tribologists and Lubrication Engineers 
Chicago, Illinois, May 14-19, 1995 



National Aeronautics and 
Space Administration 


(NASA— TM— 106585) COMP AR ISDN OF 
LIFE THEORIES FOR ROLLING-ELEMENT 
BEARINGS (NASA. Lewis Research 
Center) 28 p 


N95-267 74 


Unci as 




G3/37 0048307 


COMPARISON OF LIFE THEORIES FOR ROLLING-ELEMENT BEARINGS 

Erwin V. Zaretsky* *, Joseph V. Poplawski', 
and Steven M. Peters^ 


ABSTRACT 

Nearly five decades have passed since G. Lundberg and A. Palmgren published their 
life theory in 1947 and 1952 and it was adopted as an ANSI/ABMA and ISO standard in 
1950 and 1953. Subsequently, many variations and deviations from their life theory have been 
proposed, the most recent being that of E. Ioannides and T.A. Hams in 1985. This paper 
presents a critical analysis comparing the results of different life theories and discussing their 
implications in the design and analysis of rolling-element bearings. Variations in the stress- 
life relation and in the critical stress related to bearing life are discussed using stress fields 
obtained from three-dimensional, finite-element analysis of a ball in a nonconforming race 
under varying load. The results showed that for a ninth power stress-life exponent the 
Lundberg-Palmgren theory best predicts life as exhibited by most air-melted bearing steels. 
For a 12th power relation reflected by modem bearing steels, a Zaretsky-modified Weibull 
equation is superior. The assumption of a fatigue-limiting stress distorts the stress-life 
exponent and overpredicts life. 


SYMBOLS 

A material constant 

a semimajor axis of contact ellipse, m (in.) 

b semiminor axis of contact ellipse, m (in.) 

C D dynamic load capacity, N (lbf) 

c critical shear stress-life exponent 

D ball diameter, m (in.) 

e Weibull slope 
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’Member STLE, J.V. Poplawski & Assoc., Bethlehem, PA. 

*J.V. Poplawski & Assoc., Bethlehem, PA. 
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exponent 


life, millions of revolutions or millions of stress cycles 

characteristic life, or operating time exceeded by 36.8 percent of a group 
of bearings, millions of revolutions or millions of stress cycles 

10-Percent life, or operating time exceeded by 90 percent of a group of 
bearings, millions of revolutions or millions of stress cycles 

50-Percent life, or operating time exceeded by 50 percent of a group of 
bearings, millions of revolutions or millions of stress cycles 

length of race, m (in.) 

life, millions of stress cycles 

Hertz stress-life exponent 

equivalent bearing load, N (lbf) 

load-life exponent 

probability of survival, fraction or percent 
maximum Hertz stress. Pa (psi) 
stressed volume, cu. m (cu. in.) 
coordinates 
depth, m (in.) 

depth to critical shear stress, m (in.) 
depth to each stressed volume, m (in.) 
principal stress, Pa (psi) 

Von Mises stress, Pa (psi) 
critical shear stress. Pa (psi) 
maximum orthogonal shear stress. Pa (psi) 
fatigue-limiting shear stress. Pa (psi) 


in-plane shear stress, Pa (psi) 


t 45 


Subscripts 

i individual volume number 

ref reference state 

V elemental volume 

VR risk volume (Ref. 7) 


INTRODUCTION 

Bearing fatigue life analysis is based on the initiation or first evidence of fatigue 
spalling on either a bearing race or a rolling element. This spalling phenomenon is time or 
cycle dependent. Generally, the spall begins in the region of maximum shear stresses, which 
is located below the contact surface, and propagates into a crack network. Failures other than 
those caused by classical rolling-element fatigue are considered avoidable if the bearing is 
properly designed, handled, installed, and lubricated and is not overloaded (Ref. 1). However, 
under low elastohydrodynamic (EHD) lubricant film conditions, rolling-element fatigue can be 
surface or near-surface initiated with the spall propagating into the region of maximum shear 
stresses. 

If a number of apparently identical bearings are tested to fatigue at a specific load, 
there is a wide dispersion in life among the various bearings. For a group of 30 or more 
bearings the ratio of the longest to the shortest life may be 20 or more (Ref. 1). 

The term “basic rating life,” as used in bearing catalogs, usually means the fatigue life 
exceeded by 90 percent of the bearings or the time before which 10 percent of the bearings 
fail. This basic rating life is referred to as “L 10 life” (sometimes called B 10 life or 10-percent 
life). The 10-percent life is approximately one-seventh of the mean life, or MTBF (mean time 
between failures), for a normal life-dispersion curve (Ref. 1). 

Weibull (Ref. 2) postulated that the dispersion in fatigue life for a homogeneous group 
of rolling-element bearings can be expressed according to the 
following relation: 

In In _L = e In — 

S L p (1) 

Equation (1) relates the probability of survival S and the life L in millions of inner- 
race revolutions. When In In (1 IS) is used as the ordinate and In L as the abscissa, bearing 
fatigue data plot as a straight line. 
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Weibull (Ref. 2) further stated that the probability of survival S could be expressed as 


In I ~ x c N e V 

S (2) 

where V is the volume representation of stress concentration referred to herein as “stressed 
volume,” x is the critical shear stress, and c is a critical shear stress-life exponent where from 
Eq. (2) 


N ~ x c/e (3) 

Lundberg and Palmgren (Ref. 3) extended the theoretical work of Weibull (Ref. 2) and 
showed that the probability of survival 5 could be expressed as a power function of 
orthogonal shear stress x 0 , life N, depth to maximum orthogonal shear stress z 0 , and stressed 
volume V. That is. 








(4) 


and where 


V = 


(5) 


Then 


1 x N e al 
In - ~ — 

S hr 1 (6) 

^ O 

Lundberg and Palmgren based their life theory on subsurface-originated fatigue, where 
high EHD film conditions would exist. However, industry practice has evolved to use 
Lundberg-Palmgren theory for both high and low EHD film conditions. Hence, the theory is 
applicable to both surface- and subsurface-originated fatigue. 


Lundberg and Palmgren (Ref. 3) obtained the following additional relation: 



(7) 


where C D , the basic dynamic load capacity, is defined as the load that a bearing can carry for 
one million inner-race revolutions with a 90-percent probability of survival, P eq is the 
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equivalent bearing load, and p is the load-life exponent. Endurance tests on ball bearings have 
shown that p = 3 (Ref. 4). Most bearing manufacturers, the American National Standards 
Institute (ANSI)/ Anti-Friction Bearing Manufacturers Association (AFBMA), and the 
International Organization for Standardization (ISO) have adopted this value of p for ball 
bearings and 10/3 for roller bearings (Refs. 5 and 6). Substituting bearing geometry and Hertz 
stress for a given load and life into Eq. (6), with appropriate exponents, results in a value of 
C D for a given bearing. 

Ioannides and Harris (Ref. 7) proposed a discrete-stressed- volume approach based on 
Lundberg-Palmgren (Ref. 3) for bearing life prediction. In their analysis they also proposed 
the concept of a fatigue limit. Zaretsky (Ref. 8) proposed a generalized Weibull-based 
methodology for structural fatigue life prediction that is also based on a discrete-stressed- 
volume approach. There are, however, distinct differences in approach between the method 
proposed by Zaretsky (Ref. 8) and that of Ioannides and Harris (Ref. 7). 

Two stress-life exponents are referred to. The critical shear stress-life exponent c is 
related to the Lundberg-Palmgren parameter c/e (Eq. (3)) and applies if V and z 0 are fixed. 
The Hertz stress-life exponent n encompasses all the effects on life that arise when the Hertz 
stress is changed on a fixed geometry by varying the normal load. The contact or Hertz 
stress-life exponent n (which is different from the critical shear stress-life exponent c) is 
related to the load-life exponent p as follows: 


n 



To complicate matters, there appears to be an issue of what the value of p should be for 
purposes of analysis. Work reported by Parker and Zaretsky (Ref. 9) suggests that for air- 
melt-processed steels, such as used by Lundberg and Palmgren (Ref. 3), the Hertz stress-life 
exponent n is approximately 9. However, for the cleaner and more recent vacuum-processed 
steels, n “ 12. 

Varying the Hertz stress-life exponent n can significantly affect life predictions for 
bearings. Using n = 9 results in a more conservative estimate of bearing life than using n = 
12. Also, the ratio of the predicted lives is a function of load and is directly related to 
Cjy/P . Hence, in the normal operating load envelope, life may be underpredicted by a factor 
of 20 C when using the ANSI/ AFBMA standards or a ninth power Hertz stress-life exponent n. 
A proper stress-life exponent then becomes more than of mere academic interest, since a 
designer requires a reliable analytic tool to predict bearing life and performance. 

In view of the aforementioned discussion, the objectives of the work reported herein 
were (1) to evaluate and compare the various life theories for rolling-element bearings, 
specifically to point contacts; (2) to evaluate the effect of the critical shear stress-life 
exponent c on the Hertz stress-life exponent n and on bearing-life prediction; and (3) to 
determine what effect the presumption of a fatigue limit has on bearing life prediction. 
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HERTZ STRESS-LIFE RELATION 


The generally accepted relation between load and life in a rolling-element bearing is 
that life varies with the inverse cubic power of load for ball bearings (point contact) and with 
the inverse fourth power of load for roller bearings (line contact). 

Because stress is proportional to the cube root of load for ball bearings, the load-life 
relation for ball bearings can be resolved into a stress-life relation, where life is inversely 
proportional to the maximum Hertz stress to the ninth power. This relation has been generally 
accepted by ball bearing manufacturers and users. There is at least one exception where a 
manufacturer has indicated that the life of ball bearings varies inversely with the fourth power 
of load (or 12th power of Hertz stress). Nevertheless, the inverse cubic load-life relation has 
been included in the ANSI/AFBMA standards for ball bearings (Ref. 5). 

Later fatigue tests of ball bearings by several investigators tended to verify this inverse 
cubic relation. Styri (Ref. 10) presented data for two types of ball bearings. For one group of 
6207-size, deep-groove ball bearings under various radial loads from 3.5 to 17.3 kN (775 to 
3880 lbf), life was inversely proportional to load to the 3.3 power. Another group of 1207- 
size, double-row, self-aligning bearings was tested such that the maximum Hertz stress at the 
outer-race-ball contact varied from 4.0 to 5.6 GPa (580 to 810 ksi). Here it was found that 
life varied inversely with the ninth power of stress (or the third power of load). Cordiano et 
al. (Ref. 11) reported load-life data for 217-size, thrust-loaded ball bearings. The resultant 
Hertz stress-life exponents were 8.1, 9.6, and 12.6 for three lubricants, a water-glycol base, a 
phosphate ester base, and a phosphate ester, respectively. McKelvey and Moyer (Ref. 12) 
reported that, with four groups of AISI 4620, carburized-steel, crowned rollers (elliptical 
contact), fatigue life varied inversely with maximum Hertz stress to the eighth to ninth power. 
Maximum Hertz stress in these tests varied from 1.8 to 3.3 GPa (262 to 478 ksi). Townsend 
et al. (Ref. 13) surface fatigue tested three groups of case-carburized, consumable-electrode- 
vacuum-arc-remelted AISI 9310, 8.89-cm (3.5-in.) pitch diameter spur gears at maximum 
Hertz stresses of 1.5, 1.7, and 1.9 GPa (222, 248, and 272 ksi). The gears were run at 
10 000 rpm and 77 °C (170 °F). The lubricant was a superrefined naphthenic mineral oil with 
an additive package. The L 10 life varied inversely with stress to the 8.4 power, but the L 50 
life varied inversely with stress to the 10.2 power. The average Hertz stress-life exponent n 
was 9.3. 

Several other investigators (Refs. 14 to 22) have reported data with bench-type, 
rolling-element fatigue testers, rather than with full-scale bearings or gears. The L 10 lives as a 
function of maximum Hertz stress for these data are shown in Figure 1, with the data 
summarized in Table 1. From the table the maximum Hertz stress for these data ranged from 
3.7 to 9.0 GPa (526 to 1300 ksi). With the exception of Greenert’s work (Ref. 19) the stress- 
life exponents ranged from 8.4 to 12.4. The data are all for AISI 52100 and AISI M— 50 steels 
(except for the data reported by Barwell and Scott (Ref. 14), who do not state the type of 
steel). At least two sets of data were for air-melted material, three were for vacuum-degassed 
material, and one was for vacuum-arc-remelted steel. The other references do not state the 
melting process. The stress-life exponents from Greenert (Ref. 19), ranging from 15 to 19, are 
much higher than those from other published data. This lack of correlation is unexplained. 
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Lorosch (Ref. 23) fatigue tested three groups of vacuum-degassed, 7205B-size, AISI 
52100 inner races at maximum Hertz stresses of 2.6, 2.8, and 3.5 GPa (370, 406, 500 ksi), 
respectively. Each group consisted of 20 races, for a total of 60. (It is assumed that all 60 
races came from a single heat of material and were of the same hardness, but Lorosch 
(Ref. 23) does not state so in his paper.) After the test runs Lorosch examined cross sections 
of the races. He observed that only at the lowest stress, 2.6 GPa (370 ksi), did no measurable 
plastic deformation occur. At the two higher stresses plastic deformations of different 
magnitudes were measured. Lorosch divided the races at the two higher stresses into groups 
according to these magnitudes. Group A included the races with the smaller deformations, and 
group B, the races with the larger deformations. The resultant stress-life exponent of groups A 
and B was 12. However, when E.V. Zaretsky reconstituted the data for groups A and B into a 
single group, designated group AB, the stress-life exponent varied from 12 to 27 depending 
on the stress range over which the exponent was calculated. Lorosch (Ref. 23) did not 
calculate the stress reduction resulting from plastic deformation of the races, nor did he report 
on component hardness, hardness differential between the balls and races, and residual 
stresses induced during operation, all of which would also have affected his results. From 
these tests Lorosch (Ref. 23) concluded that “under low loads and with elastohydrodynamic 
lubrication there is no material fatigue, thus indicating that under such conditions bearing life 
is practically unlimited.” 

Zwirlein and Schlicht (Ref. 24), in a companion paper published concurrently with that 
of Lorosch (Ref. 23) and using the same 7205B-size bearing inner race data, state that 
“contact pressures less than 2.6 GPa (370 ksi) do not lead to the formation of pitting within a 
foreseeable period. This corresponds to ' true endurance strength ’ .” This observation would 
support the assumption by Ioannides and Harris (Ref. 7) of the existence of a fatigue limit 
for bearing steels.” However, this observation is not supported by rolling-element fatigue data 
in the open literature for stress levels under 2.6 GPa (370 ksi) — such as those reported in 
References 25 to 29, which exhibited classical rolling-element fatigue. In rotating machinery 
nearly all rolling-element bearings operate at a maximum Hertz stress less than 2.1 GPa 
(300 ksi). Therefore, if Lorosch (Ref. 23) and Zwirlein and Schlicht (Ref. 24) were correct, 
no bearing in rotating machinery applications would fail by classical rolling-element fatigue. 

To the authors’ knowledge there are no reported laboratory-generated, full-scale 
bearing rolling-element fatigue data at stresses significantly below 2.1 GPa (300 ksi) that 
would either establish or refute this presumption. However, Townsend et al. (Ref. 30) 
reported rolling-element fatigue tests for through-hardened spur gears made from vacuum- 
induction-melted, consumable-electrode- vacuum-arc-remelted (VEM-VAR) AISI M-50 steel 
run at a maximum Hertz stress of 1.7 GPa (248 ksi). The AISI M-50 gears ran for at least 
five times the life of identical AISI 9310 gears and failed both by classical rolling-element 
fatigue and tooth bending fatigue. What is significant about these tests is that classical rolling- 
element fatigue is reported for a bearing steel (M-50) at stresses as low as 1.7 GPa (248 ksi), 
the stress range in which most bearings operate. If there exists a fatigue limit for bearing 
steels, it certainly would have to be less than a maximum Hertz stress of 1.7 GPa (248 ksi). 

The explanation for the trend in the Lorosch (Ref. 23) and Zwirlein and Schlicht 
(Ref. 24) data is the inducement of compressive residual stresses in the AISI 52100 steel 
caused by the transformation of retained austenite into martensite during rolling-element 
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cycling. These compressive residual stresses were reported by Zwirlein and Schlicht (Ref. 24). 
Compressive residual stresses reduce the effective magnitude of the maximum shear stresses 
caused by Hertzian loading. This lower stress results in longer bearing life and deviation from 
the Hertz stress-life exponent n of 9 or 12 to a significantly higher value. Lorosch (Ref. 23) 
and Zwirlein and Schlicht (Ref. 24) extrapolated their data, leading them to conclude the 
existence of a fatigue-limiting stress rather than that induced compressive residual stresses had 
increased bearing life. 

Lorosch (Ref. 23) performed another series of rolling-element fatigue experiments with 
20 inner races, designated group C, from the same heat as groups A and B. He used a 
Rockwell hardness tester to make 0.1-mm-diameter indentations at four evenly spaced 
locations around the center of the race circumference. He divided group C into two 10- 
bearing sets, which were tested at 2.6 and 2.8 GPa (370 and 406 ksi), respectively. The Ljq 
lives for group C were significantly reduced from the group AB lives. These tests exhibited a 
9.5 stress-life exponent. Because these indentations are analogous to wear debris denting 
during bearing operation, the results suggest that surface damage is another factor affecting 
the stress-life exponent. 


BALL-RACE MODEL 

To compare life theories and their variants, we selected a simple ball-race geometry 
model for evaluation. The model assumes a plurality of normally loaded 12.7-mm (0.5-in.) 
diameter balls running in a linear, nonconforming grooved race having a groove radius of 
6.6 mm (0.26 in.) and a length /. A schematic of the ball-race model is shown in Figure 2. 
The loads, stresses, and dimensions used for the ball-race model of Figure 2 are summarized 
in Table 2. The three-dimensional, finite-element analysis (FEA), geometry-for-the-gutter 
(Fig. 2) model used in the studies is shown in Figure 3. The model geometry takes advantage 
of the symmetric nature of the Hertzian contact for the case of no significant surface shear 
stresses. This model was expanded to encompass the full contact ellipse for our additional 
studies, which included different surface shear patterns and their influence on the related 
stress fields. However, those results are not reported as a part of this study. A full-race FEA 
geometry was also constructed. The full-race model will be used for studying the interaction 
of the centrifugal, axial clamping and press-fit stresses on the three-dimensional contact stress 
fields and their influence on bearing fatigue life. 

The quarter section of the contact ellipse face was divided into approximately 162 
elements. Element size ranged from 0.0991 mm x 0.0330 mm (0.0039 in. x 0.0013 in.) to 
0.1278 mm x 0.03175 mm (0.00503 in. x 0.00125 in.) depending on the Hertzian stress level. 
Element thickness in the depth direction was 0.0254 mm (0.0010 in.) until a depth zJb of 
about 1.0. Beyond that depth the element thickness was gradually increased. A typical model 
contained approximately 3500 to 5900 solid isoparametric elements depending on the Hertz 
stress. The model for 2.4-GPa (350-ksi) Hertz stress had approximately 5800 elements and 
7000 nodes, giving about 18 000 degrees of freedom after applying constraint boundary 
conditions. The analysis was performed on a 386/33-MHz personal computer with the 
COSMOS/M commercially available FEA software. 
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We checked the FEA model results against calculated values by using classical Hertz 
contact stress theory. Figure 4 shows a comparison of the subsurface stresses at the center of 
the contact ellipse for 1.4-GPa (200-ksi) Hertzian stress. The FEA-predicted principal normal 
stresses and the in-plane shear stress x 45 agreed within 3 percent of theory over the Hertz 
stress range studied. 

Figure 5 shows the three-dimensional Von Mises effective or equivalent stress field. A 
maximum stress of about 0.84 GPa (122 ksi) occurred 0.114 mm (0.005 in.) below the 
surface for a 1.4-GPa (200-ksi) Hertz stress. Figure 6 shows that the in-plane shear stress X 45 
on a plane cut through the center of the contact ellipse perpendicular to the rolling direction. 
The maximum shear stress of about 0.45 GPa (66 ksi) also occurred 0.114 mm (0.005 in.) 
below the surface. Finally, Figure 7 shows the maximum orthogonal shear stress field acting 
in the y-z plane. A mirror image of \ yz also occurred at the forward half of the vertical 
symmetry axis but had the opposite sign (not shown). The largest orthogonal shear stress was 
0.32 GPa (46 ksi). It was the double amplitude of x yz on which Lundberg and Palmgren based 
their selection of that stress for their analysis. 

Figures 5 to 7 were included to show the three stress distributions that have been 
discussed over the years as being the critical stress in determining bearing fatigue life. These 
three stresses are examined as the “stress of choice” within this paper. 

The results of the FEA runs at the three Hertz stresses were saved as data bases to be 
used in evaluating the life theories examined in this paper. For purposes of analysis, only the 
life of the race will be considered at each load condition. The normal loads per ball are 623, 
1617, and 3333 N (140, 364, and 749 lbf) producing maximum Hertz stresses of 1.4, 1.9, and 
2.4 GPa (200, 275, and 350 ksi), respectively. The L 10 life at 2.4 GPa (350 ksi) is 
normalized and assumed to be 1. The normalized fives at each stress load are summarized in 
Table 2 for stress-fife exponents n of 9 and 12 and shown in Figure 8, where 


(9) 


RESULTS AND DISCUSSION 


Weibull Equation 

The Lundberg-Palmgren theory (Ref. 3) is a variant of the Weibull theory of failure 
(Ref. 2). The genius of Lundberg and Palmgren was to realize its applicability to rolling- 
element bearings and to fit both bearing geometry and fife data to the Weibull parameters. 
However, Lundberg and Palmgren do deviate from the Weibull equation by introducing the 
depth to the maximum critical shear stress z 0 shown in Eqs. (4) and (5), which they assumed 
to be the orthogonal shear stress. The rationale for introducing z 0 was, according to them, 
because “it is possible that changes in condition of the material also depend on the depth 
below the surface of the volume element considered. This assumption is necessary for the 
treatment of fatigue in rolling bearings.” They go on to state, “Weibull’ s theory, however, is 
built on the assumption that the first crack leads to a break. The frequent examples of cracks 
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which do not reach the surface seem to show that this assumption is not valid for fatigue in 
rolling bearings. Therefore, when setting up the life formulas consideration ought to be given 
to the fact that the probability of a fatigue break occurring must be taken to be dependent on 
the depth z 0 at which the most dangerous stress occurs” (Ref. 3). 

Lundberg and Palmgren did not consider that in high-cycle fatigue crack initiation is 
probably over 95 percent of the life cycle and crack propagation is probably less than 5 
percent of the life cycle. Accordingly, the life formula without z 0 , as proposed by Weibull 
(Eq. (2)), may be applicable to bearing life prediction. Equation (2) can be written as follows: 


LlO - A 




1 r ' c 

- ( i 




max 


( 10 ) 


The factor A is a material constant. From Hertz theory V and x can be expressed as a function 
of S max where 


x 


max 


( 11 ) 


V - 


Substituting these values into Eq. (10) and solving for the exponents give 


( 12 ) 


Experience has shown that the Weibull slope e for most bearing fatigue data varies from 1 to 
2. However, it has been accepted practice since Lundberg-Palmgren (Ref. 3) to assume this 
value to be 1.1. Experience has also shown the value of the Lundberg-Palmgren parameter c/e 
to be 9. As a result, n would equal 10.81. 

Assuming n of 9 and 12, c/e would equal 7.18 and 10.18, respectively. Whereas c/e = 
7.18 can be considered too low, c/e = 10.18, although slightly higher than 9, may be 
acceptable. Normalized values of the material constant A for a Hertz stress-life exponent n of 
10.81 are given in Table 3. The calculated relative lives for each of the stresses are given in 
Table 4. 

Lundberg-Palmgren Equation 

The Lundberg-Palmgren theory has been discussed in some detail. A reading of 
Reference 3, however, does not make intuitively obvious their approach to bearing life 
prediction and the determination of their variables and exponents. Equation (4) can be 
rewritten as follows: 
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max 


( 14 ) 


As for V and x 0 from Eqs. (11) and (12), z 0 can be expressed as a function of S max , where 


Substituting these values into Eq. (14) and solving for the exponents give 

-2 - c + h _ 

7 ' ( 16 ) 

If it is assumed that g = 1.1 and n = 9, there are two unknowns to be solved for. A second 
equation is therefore necessary. 

Lundberg and Palmgren (Ref. 3) showed rolling-element fatigue data (L 50 ) from which 
a relation between ball diameter D and life was empirically determined, where 

L - D 5a 

(17) 

From Hertz theory a relation for V, % Q , z 0 , and D may be obtained as follows: 


max 


(15) 


-2 


(18) 


V ~ azj ~ abl ~ D 3 D ~ D 


2 


5 

3 


(19) 



( 20 ) 


Substituting Eqs. (17) to (20) into Eq. (14) and solving for the exponents give 


-5 + 2c + h 


( 21 ) 
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Combining Eqs. (21) and (16) and substituting e = 1.1 and n = 9 yield h = 2.33, or hJe - 2.1 
and c/e = 9.3. These values, summarized in Table 3, are those of the Lundberg-Palmgren 
equation (Ref. 3) and the ANSI/AFBMA standards (Ref. 5). 

What Lundberg and Palmgren did not consider was the statistical variability of their 
data related to Eq. (16). From Hertz theory 


'max 


- D 


-213 


( 22 ) 


Substituting Eq. (22) together with Eqs. (18) to (20) into Eq. (14) and solving for the 
exponents give 


-5 + 2c + h _ In 

3e " T (23) 

Letting e = 1.1 and n = 9 yields h = 3, or hie = 2.7 and c/e = 9.9. Although c/e = 9.9 is 
slightly high, it is acceptable. If it is assumed that n = 12, Eqs. (16) and (23) yield hie = 2.7 
and c/e = 12.91. From experience this c/e is too high and not realistic. 

Ioannides-Harris Equation 

Ioannides and Harris (Ref. 7), using Weibull (Ref. 2) and Lundberg and Palmgren 
(Ref. 3), introduced a fatigue-limiting stress and integrated the computed life of elemental 
stress volumes to predict bearing life. Their equation, which relates to each elemental volume, 
is as follows: 


1 1 

In - ~N Jvt? t— 


dV 


(24) 


The equation is similar to that of Lundberg and Palmgren (Eq. (4)) except for the introduction 
of a fatigue-limiting stress X u . Equation (24) can be rewritten to represent each individual 
stressed volume as follows: 


L v - A 


1 

J~V , 

J 


1 ( 
e 


T - X 

v 0 “/ 


(£y) 


(25) 


The three-dimensional FEA model data base previously discussed contains the elemental 
volumes and associated stresses. An elemental reference volume is required to be selected as 
well as a reference stress and reference depth. These reference points can be taken as the 
volume comprising the maximum critical shear stress and the depth to this volume. Equation 
(25) can be rewritten as follows: 


12 


( 26 ) 


V 


1 

1 

e 

(^0 ~ ^1 ) ref 

C 

e 

Z/ 

h 

e 

J- ref_ 


[ v i \ 


1 

1 

1 


^ref 




1 


v 



/ 


(27) 


Equation (26) is identical to the Lundberg-Palmgren Eq. (14) except for the fatigue-limiting 
stress x u . Values of e = 1.1, c/e = 9.3, and hie = 2.1 for the ball-race model of Figure 2 as 
well as a reference volume and shear stress T45 of 104x10 ^ m^ (6.33x10 ^ in.^) and 0.77 
GPa (111 ksi), respectively, were assumed for the Ioannides-Harris Eq. (26). Values of x u = 
0.21 GPa (30 ksi) for x 0 and x 45 and 0.35 GPa (50 ksi) for were also assumed. The 
relative life calculations are summarized in Table 4 and compared with the previous 
calculations for the Weibull Eq. (10) and the Lundberg-Palmgren Eq. (14). As might be 
expected, the relative life values at 1.4 GPa (200 ksi) are extremely high with this approach, 
in contrast to what would be expected. If the apparent values of the Hertz stress-life exponent 
n were 19.6 when the orthogonal shear stress was used, 16.8 when the maximum shear stress 
was used, and 15.1 when the Von Mises stress was used, only two sets of data in the 
literature (Refs. 19, 23, and 24) would support this analysis. 


Zaretsky Equation 

Both the Weibull and Lundberg-Palmgren equations relate the critical shear stress-life 
exponent c to the Weibull slope e. The Lundberg-Palmgren parameter de thus becomes, in 
essence, the effective critical shear stress-life exponent, implying that the critical shear stress- 
life exponent depends on bearing life scatter or dispersion. A search of the literature for a 
wide variety of materials and for non-rolling-element fatigue reveals that most stress-life 
exponents vary from 6 to 12. The exponent appears to be independent of scatter or dispersion 
in the data. Hence, Zaretsky (Ref. 31) has rewritten the Weibull equation to reflect that 
observation by making the exponent c independent of the Weibull slope e, where 


In A - N e V 

S 


(28) 


Equation (28) differs from the Weibull Eq. (2) and the Lundberg-Palmgren Eq. (4) only in the 
exponent of the stress x. Only for e = 1 do the two reflect the same power relation. Equation 
(28) differs from Lundberg-Palmgren further in that it does not recognize the effect of stress 
depth z 0 . In addition, the Zaretsky model differs from Lundberg-Palmgren conceptually 
because it makes survival probabilities dependent on local stress pointwise in the stressed 
volume as does Ioannides-Hams (Ref. 7), whereas Lundberg-Palmgren used only one critical 
stress. Equation (28) can be written like Eq. (10), 


L10 - A 




fi y 

A, 


(29) 
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Equation (29) has a life-stressed-volume exponent dependent on data dispersion or cumulative 
distribution and a critical stress exponent independent of data dispersion or Weibull slope. 

Zaretsky, like Ioannides and Harris (Ref. 7), integrates the complete life of elemental 
stressed volumes as follows: 


In 1 ~ N e L dV 

S (30) 

As with Ioannides and Harris (Ref. 7) an elemental reference volume is required. Equation 
(29) is rewritten as follows: 


“ - 


r _ _ i 

Vie 

~ ~ 

L i 


*ref 


X ref 

_V] 




1 


The lives of the individual stressed elements are summarized in accordance with Eq. (27). 

Using the same finite-element analysis as used for the Ioannides-Harris Eq. (26) and 
assuming that c = 9 and e = 1.1, the relative lives at each stress level were calculated and are 
summarized in Table 4 on the basis of maximum shear stress x 45 and the Von Mises stress 
<7^. The apparent Hertz stress-life exponent n was 10.8 for the maximum shear stress and 
11.1 for the Von Mises shear stress. Both values coincide with experimental values in the 
literature (Table 1). However, the same stress-life exponent is obtained if the Weibull Eq. 

(10) is used without a finite-element analysis. 

Although Zaretsky (Ref. 31) does not propose a fatigue-limiting stress, he does not 
exclude that concept either. However, his approach is entirely different from that of Ioannides 
and Harris (Ref. 7). For critical stresses less than the fatigue-limiting stress, the life for the 
elemental stressed volume is assumed to be infinite. With this assumption and for the same 
fatigue-limiting stresses used previously with the Ioannides-Harris Eq. (26), the lives at the 
three stress levels were recalculated. The resultant fives were no different from those shown 
under the Zaretsky equation column in Table 4. However, if (x i - x u ) and (x ref - x u ) are 
substituted for x, and x ref in Eq. (31), the values of the stress-life exponent n approach those 
of Ioannides and Harris (Ref. 7) shown in Table 4. 

MATERIAL CONSTANT 

The material constant A in each life equation allows the equation to be matched with 
existing data, as Lundberg and Palmgren (Ref. 3) did in their original work. Because there are 
no data for the ball-race model of Figure 2, it is possible to normalize the material constant 
for a normalized L 10 fife of 1 at a maximum Hertz stress S may of 2.4 GPa (350 ksi). These 
normalized values are summarized and compared in Table 3 for each fife analysis. The units 
for A will differ depending on what equation it is related to. The race length l is factored out 
of the stressed volume and is shown as part of the material constant. For finite-element 
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analyses the length of the unit volume in the rolling direction has been factored into the 
material constant. 


CONCLUDING COMMENTS 

The comparisons presented herein are perhaps the first time a qualitative comparison 
has been made between the various theories and approaches to predicting bearing life. No 
attempt was made to factor into the life equation other variables covered by STLE Life 
Factors for Rolling Bearings (Ref. 32). However, it is noteworthy to consider the effect of 
residual stresses and hoop stresses on rolling-element fatigue life. The Lundberg-Palmgren 
equation uses orthogonal shear stress as the critical shear stress. This stress would not be 
significantly affected either by hoop or residual stresses below the contact surface. The two 
dominant stresses that must be considered are the shear stress x 45 and the Von Mises stress 
O . Our calculations strongly suggest that the Von Mises stress may more accurately reflect 
the magnitude of change in life with hoop stress and that both the Von Mises and T 45 shear 
stresses may proportionately reflect the effect of residual stress. 

We provided a finite-element method that permits the practical implementation, with 
accessible computer resources, of a life model which integrates survival probabilities over the 
stressed volume for several suggested expressions for the survival probability of a volume 
element, such as those of Lundberg-Palmgren, Ioannides-Harris, or Zaretsky. Using this 
method it should now be possible to more accurately assess the effect of surface traction, 
residual stress, and hoop stress on bearing life and reliability. 

In our opinion there may be strong justification that this is the time to reexamine the 
Lundberg-Palmgren equation. The equation has served the manufacturer, designer, and user of 
bearings well. However, even with STLE life factors the life predictions are conservative. We 
believe that the approach of Zaretsky (Ref. 31), based on a modified Weibull equation, best 
reflects the modem data base for bearing steels. However, a significant effort would be 
required to change existing bearing design computer codes and the ANSI/AFBMA and ISO 
standards in order to make any changes practical to the user community. 

SUMMARY OF RESULTS 

A critical analysis was presented comparing the results of different life theories and 
discussing their implications in the design and analysis of rolling-element bearings. The 
following results were obtained: 

1 , For an inverse ninth power relation between life and maximum Hertz stress, the 
Lundberg-Palmgren theory qualitatively predicts life best. However, for an inverse 12th power 
relation between life and maximum Hertz stress, a Zaretsky-modified Weibull equation is 
superior. 

2. The Ioannides-Harris-modified, Lundberg-Palmgren equation based on a “fatigue- 
limiting stress” significantly overpredicts life. Analysis implies an inverse life— maximum 
Hertz stress exponent as high as 19.6 for a fatigue-limiting stress of 0.21 GPa (30 ksi). 
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3. Using a fatigue-limiting stress coupled with a finite-element analysis while 
assuming that those elements being stressed below the fatigue limit have infinite life does not 
affect life prediction. However, subtracting a fatigue-limiting stress from each elemental stress 
distorts the life-maximum Hertz stress relation, producing unreasonably high predictions of 
life. 
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TABLE 1.— PUBLISHED BENCH-TYPE RIG ROLLING-ELEMENT FATIGUE DATA RELATED TO STRESS-LIFE EFFECTS 

[From Parker and Zaretsky (Ref. 9).] 


Number 
from 
Fig. 1 

Reference 

Material 

Lubricant 

Load range, 
kg 

Maximum Hertz 
stress range, 
GPa (ksi) 

Load- 

life 

exponent 

Stress- 

life 

exponent 

Test type 


Type 

Melting 

process 



1 

Barwell and 
Scott (Ref. 14) 


(a) 

Mineral oil 

200 to 600 

b 4.5 to 6.6 ( b 660 to 950) 

2.8 

8.4 

Four ball 

2 

Butler and Carter 
(Ref. 15) 

AISI 52100 

Air melt 

SAE 10 
mineral oil 


4.1 to 5.2 (600 to 750) 

3.5 

10.4 

Spin rig 

3 

Butler and Carter 
(Ref. 15) 

MV-1 

(AISI M-50) 

Air melt 

SAE 10 
mineral oil 


4.1 to 5.2 (600 to 750) 

3.2 

9.7 

Spin rig 

4 

Baughman 
(Ref. 16) 

MV-1 

(AISI M-50) 

(a) 

MII^k-7808 


4.4 to 5.4 (640 to 777) 

3.2 

9.7 

Rolling- 
contact rig 

5 

Scott (Ref. 17) 

EN-31(AISI 

52100) 

(a) 

Diester 

400 to 600 

b 5.7 to 6.6 (830 to 950) 

3.6 

10.8 

Four ball 

6 

Ulemi and 
Okomoto 
(Ref. 18) 

AISI 52100 

(a) 

#60 Spindle oil 


3.7 to 5.1 (526 to 730) 

00 

c4 

o 

c 8.5 

Crowned 

disk 

7 

Greenert 
(Ref. 19) 

AISI 52100 

(a) 

Navy 2190 TEP 
lubricating oil 


5.8 to 6.1 (840 to 890) 

c 5 to 6.3 

c 15 to 19 

Toroids 

8 

Valori et al. 
(Ref. 20) 

AISI 52100 

Vacuum-arc 

remelted 

Mineral oil 


4.2 to 5.5 (610 to 800) 

d 4.1 

d 12.4 

Four ball 

9 

Schatzberg and 
Felsen (Ref. 21) 

AISI 52100 

Vacuum 

degassed 

Squalene 


6.4 to 9.0 (925 to 1300) 

3.8 

11.5 

Four ball 

10 

Schatzberg and 
Felsen (Ref. 21) 

AISI 52100 

Vacuum 

degassed 

Squalene + 
100 ppm H 2 0 


6.4 to 9.0 (925 to 1300) 

3.9 

11.7 

Four ball 

11 

Parker and 
Zaretsky (Ref. 9) 

AISI 52100 

Vacuum 

degassed 

Paraffinic 
mineral oil 


4.5 to 6.0 (650 to 857) 

4 

12 

Five ball 


a Melting process not reported. Estimated, not reported by authors of reference. 

Approximate stress range, not reported by authors of reference. Reanalyzed by Zaretsky and Parker (Ref. 22). 


TABLE 2 — SUMMARY OF STRESSES, HERTZIAN DIMENSIONS, AND RELATIVE LIVES FOR 
BALL-RACE MODEL OF FIGURE 2 
[Ball diameter, 12.7 mm (0.5 in.); race-groove radius, 6.6 mm (0.026 in.).] 


Parameter 

Normal ball load, N (lbf) 


623 (140) 

1617 (364) 

3333 (749) 

Maximum Hertz stress, GPa (ksi) 

1.4 (200) 

1.9 (275) 

2.4 (350) 

Hertzian contact ellipse, mm (in.): 




Semimajor axis, a 

1.3 (0.05) 

1.8 (0.07) 

2.3 (0.09) 

Semiminor axis, b 

0.16 (0.006) 

0.22 (0.009) 

0.28 (0.011) 

Maximum shear stress values and depths: 




V GPa (ksi) 

0.343 (49.8) 

0.472 (68.5) 

0.601 (87.2) 

z Q , mm (in.) 

0.079 (0.003) 

0.107 (0.004) 

0.137 (0.005) 

W GPa (ksi) 

0.459 (66.6) 

0.644 (93.4) 

0.814(118.1) 

z Q , mm (in.) 

0.114(0.005) 

0.140 (0.006) 

0.191 (0.008) 

Ovnr GPa (ksi) 

0.821 (119.1) 

1.154 (167.4) 

1.465 (212.4) 

z 0 , mm (in.) 

0.114(0.005) 

0.140 (0.006) 

0.191 (0.008) 

Relative lives for values if n from Eq. (9): 




n = 9 

153.94 

8.76 

1.0 

n = 12 

825.01 

18.06 

1.0 


TABLE 3.— ASSUMED AND CALCULATED VALUES FOR LIFE ANALYSIS OF TABLE 4 


Parameter 

Weibull 

Lundberg- 

Ioannides-Harris Eq. (25) 

Zaretsky 


Eq. (10) 

Palmgren 




Eq. (29) 



Eq. (14) 






Critical shear stress. 

X max 

% 


*45 

°vm 

X 45 

°vm 

GPa (ksi) 








Weibull slope, e 

1.1 

u 

1.1 

1.1 

1.1 

1.1 

1.1 

Lundberg-Palmgren 








parameters: 






9.0 b 


de 

9.0 

9.3* 

9.3 

9.3 

9.3 

9.0 b 

hie 


2.1* 

2.1 

2.1 

2.1 



Shear fatigue- 
limiting stress, z u 



0.21 (30) 

0.21 (30) 

0.35 (50) 







Normalized material 

2.6x1 0 74 

8.1xl0 80 

0.58x1 0 39 

4.5x1 0 39 

2.5x10 s9 

1.5X10 38 

1.3X10 36 

constant, 4 A// 0909 , 
N-m (lbf- in.) 

(5.8X10 42 ) 

(1.6X10 46 ) 

(11 542) 

(90 356) 

(49 358) 

(34 056) 

(28 357) 

Stress-life exponent. 

10.8° 

9.0 

19.6 C 

16.8 C 

15.1 c 

10.8° 

11.1® 

n (Bq. (9)) 









Calculated. 

kValue of exponent c where Lundberg-Palmgren parameter de becomes cele. 
Calculated (refer to Table 4). 
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TABLE 4.— COMPARISON OF RELATIVE LIFE FROM FOUR LIFE EQUATIONS 


Maximum 
Hertz stress, 
GPa (ksi) 

Relative lives 
for values of 
n from Eq. (9) 

Relative theoretical lives 

n = 9 

n = 12 

Wei bull 
Eq. (10)“ 

Lundberg- 
Palmgren 
Eq. (14) b 

Ioannides-Harris Eq. (25) 

Zaretsky Eq. (29) 

Based on 
*o 

Based on 
*45 

Based on 

°vm 

Based 
°n *45 

Based 

ono vm 

1.4 (200) 

153.9 

825.0 

423.9 

154.7 

59 477.8 

12 112.0 

4643.0 

431.4 

495.6 

1.9 (275) 

8.8 

18.1 

13.6 

7.1 

76.1 

149.7 

85.2 

11.7 

12.9 

2.4 (350) 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Calculated 

9.0 

12.0 

c 10.8 

^.0 

c 19.6 

c 16.8 

c 15.1 

c 10.8 

c ll.l 

value of n 











a Based on x^. 
hBased on x Q . 

c Determined between 1.4 and 2.4 GPa (200 and 350 ksi). 



Life, millions of stress cycles 

Figure 1 . — Summary of published stress-life relation data for Hertzian 
contacts failing from rolling-element fatigue. (Refer to Table 1 .) Ref. 9. 
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Normal load 



Figure 3. — Finite-element quarter symmetry ball-race 
model. 
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Normalized stress, t a JS or o/S 



Normalized depth below groove surface, Z/b 


Figure 4. — Comparison of finite-element stress analysis (FEA) with closed-form solution for 
principal stresses and t 45 shear stress as a function of depth below groove surface (see 
Fig. 2). Ball diameter, 12.7 mm (0.5 in.); groove conformity, 52 percent; normal load, 623 N 
(140 Ibf); maximum Hertz stress, 1 .4 GPa (200 ksi); semimajor contact axis, a, 1 .3 mm 
(0.05 in.); semiminor contact axis, 6, 0.16 mm (0.006 in.). 
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Von Mises stress, 

rr vnT 
GPa (ksi) 

0.841 (122) 

0.731 (106) 

0.629 (91 .2) 

0.419(60.8) 

0.314(45.6) 

0.210(30.5) 

0.001 (0.133) — 



Figure 5. — Von Mises stress field in race groove for Hertzian contact of ball on groove (see 
Fig. 2). Ball diameter, 12.7 mm (0.5 in.); groove conformity, 52 percent; normal load, 623 N 
(140 Ibf); maximum Hertz stress, 1 .4 GPa (200 ksi); semimajor contact axis, a, 1 .3 mm 
(0.05 in.); semiminor contact axis, b , 0.16 mm (0.006 in.). 
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In-plane 
shear stress, 


r 45’ 


GPa 

(ksi) 

0.452 

(65.5) 

0.421 

(61.1) 

0.392 

(56.8) 

0.331 

(48.1) 

0.301 

(43.7) 

0.271 

(39.3) 

0.241 

(35.0) 

0.211 

(30.6) 

0.181 

(26.3) 

0.151 

(21.9) 

0.121 

(17.5) 

0.053 

(7.8) 

0.091 

(13.2) 



z 

A 

y 


■►x 


Figure 6. — In-plane shear stress in race groove for Hertzian contact of ball on groove (see Fig. 2). 

Ball diameter, 12.7 mm (0.5 in.); groove conformity, 52 percent; normal load, 623 N (140 Ibf); maximum 
Hertz stress, 1 .4 GPa (200 ksi); semimajor contact axis, a, 1 .3 mm (0.05 in.); semiminor contact axis, b, 
0.16 mm (0.006 in.). 
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Orthogonal 
shear stress, 


GPa (ksi) 



Figure 7. — Orthogonal shear stress in race groove for Hertzian contact of ball on groove (see Fig. 2). Ball diameter, 

12.7 mm (0.5 in.); groove conformity, 52 percent; normal load, 623 N (140 Ibf); maximum Hertz stress, 1.4 GPa (200 ksi); 
semimajor contact axis, a, 1.3 mm (0.05 in.); semiminor contact axis, b, 0.16 mm (0.006 in.). 
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Figure 8. — Effect of stress-life exponent and maximum Hertz stress on rolling-element 
fatigue life of nonconforming race groove. 
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